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ABSTRACT 

The collisionless dissipation of anisotropic Alfvenic turbulence is a promising candidate to 
solve the solar wind heating problem. Extensive studies examined the kinetic properties of Alfven 
waves in simple Maxwellian or bi-Maxwellian plasmas. However, the observed electron velocity 
distribution functions in the solar wind are more complex. In this study, we analyze the properties 
of kinetic Alfven waves in a plasma with two drifting electron populations. We numerically solve 
the linearized Maxwell-Vlasov equations and find that the damping rate and the proton-electron 
energy partition for kinetic Alfven waves are significantly modified in such plasmas, compared 
to plasmas without electron drifts. We suggest that electron drift is an important factor to take 
into account when considering the dissipation of Alfvenic turbulence in the solar wind or other 
d ~ 1 astrophysical plasmas. 


1. Introduction 


One major question in space physics is how 
collisionless plasmas in the extended corona and 
in the solar wind are heated. Observations in the 
solar wind reveal that the fluctuations of the elec¬ 
tric and magnetic fields show a turbulent spectrum 
similar to the power-law spectrum of fluid turbu¬ 
lence as described by Kolmogorov (ITu &: Marsch 


1995: Bale et al 


Alexandrova et al 


20051: 

201.C 


Bruno &: Carbone 2013t 


The turbulence in 


the s olar wind shows main ly Alfvenic polariza¬ 
tion ( Belcher fc Davij 197ll) and becomes more 
anisotr opic during the cascade to h i gher wavenum¬ 
bers fc (iGoldreich fc Sridh^ll995l : iHorburv et al 


20081 : IChen et ^ 12010 ). resulting in fc_L 
at short wavelengths. At the proton scale, the 
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turbulence is thought to tran sition into a kinetic 


Alfve n wave (KAW) cascade ([Schekochihinet_^ 
2009; Sahraoui et al. 2010: Howes et al. 2011 


Salem et al.l2012t Boldvrev fc Perej2012t Chen et al. 


20131) . w hich gradually dissipates energy to th e 
particles ( Leamon et al.lll9^ iHowes et al.ll2011 ). 


Despite the nonlinear nature of turbulence, 
linear theory has been used in studies of solar 


wind turbulence (see, for instance, I Howes et al 


(120141) for review). These studies analyze turbu¬ 
lent fluctuations by means of the linear propa¬ 
gation and damping characteristics of KAWs to 
help to understand the dissipation of the turbu¬ 
lence and the heating of the solar wind. The 
vast majority of these studies rely on represent¬ 
ing the velocity distribution functions (VDFs) 

by an isotropic or _ bi- Maxwellian bac k groun d 

(e.g., Quataert (1998 ); Leamon et al.l ( 1999f ): 
ICranmer fc van BallegooiiCT ( 20031 ) 1. In this 
work, we investigate how more realistic VDFs 
modify the damping of KAWs. In the solar 
wind, both ions and electrons can be modeled 
by several populations drifting with respect to 
each other. Recent studies show that the differ- 
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ential flow between ions affects ion heatin g by 
both ion-cyclotron wa ves ( Kas^e^^et^jiy 20^^ and 
low-frequency KAWs ( Chandran et al.ll2013h . 

The solar-wind electron VDF can be modeled 
as a superposition of three electron populations: a 
cool and dense “core” , a hot and less dense “halo” 
( Feldman et al.| 1975 ) and a bea m-like one-sided 
“strahl” ( Maksimovic et aT1 12005 b In the proton 
frame, core (halo and strahl) electrons drift sun¬ 
ward (anti-sunward) along the background mag¬ 
netic field. Empirically, the solar wind fulfills 
quasi-neutrality and the zero-current co ndition 


( Feldman et al. 19751: IPuhioa et al. 20l4) . The 


uncertainty in the fitting parameters for electron 
core properties is smaller than the uncertainty in 
the fitting parameters for halo and strahl proper¬ 
ties. A drifting bi-Maxwellian population mod¬ 
els core electrons very well. In particular, the 
core electron bulk drift speed is usually compara¬ 
ble to or larger than the Alfven speed and shows 
a clear statistical d ependence on collisional age 
( Pulupa et al.l 12014 b In the extrapolated asymp¬ 
totic limit of no collisions, the core drift is as large 
as three to four times the Alfven speed. 


In this letter, we examine how the drift between 
electron core and halo affects the linear damping 
of KAWs. We model the solar wind electrons 
by a superposition of a Maxwellian core and a 
Maxwellian halo drifting against each other in the 
proton frame. Choosing a set of parameters typi¬ 
cal for the solar wind, we solve the full hot-plasma 
dispersion relation of KAWs in the framework of 
Vlasov-Maxwell theory and identify the contribu¬ 
tion to wave damping from each plasma compo¬ 
nent. We show that, in the linear approximation, 
electron drifts can lead to significant variations in 
both the KAW damping rates and the relative en¬ 
ergy transfer from the waves to protons and elec¬ 
trons, the latter of which is another important un¬ 
solved problem in turbulent plasma heating. 


2. Theory and method 

Motivated by observed solar wind electron 
VDFs, we consider a collisionless, homogeneous, 
and warm proton-electron plasma in a constant 
uniform background magnetic field. The proton 
temperature and density are denoted as Tp and 
Up; the electron distribution consists of two shifted 
Maxwellian populations, a cool (Tc) “core” and a 


hot (T/j) “halo.” As an initial step to study the 
effects of core electron drifts on resonant damping, 
we use isotropic temperatures for all plasma com¬ 
ponents. Both electron populations drift along 
the background magnetic field Bq = BqZ with 
Vc = Vci and v/j = vhi in the proton frame. The 
halo electron drift adjusts to the core electron 
drift to guarantee the absence of net currents in 
the plasma: Vh = —VcTic/nh, where tic and nh 
denote core and halo electron densities. Figured] 
shows the plasma described above schematically. 

Under the assumption of strong anisotropy in 
kinetic scale fluctuations in the solar wind, we con¬ 
sider KAWs with wavevectors at a large angle with 
respect to the magnetic field. We use constant val¬ 
ues for the angle between k (wavenumber vector) 
and Bq (89° or 91°), plasma beta, the electron- 
core-to-proton temperature ratio and the core-to- 
halo density ratio, so that the only free parame¬ 
ters are k, Vc and sign(k • Bq). In particular, we 
allow Bq to be either parallel or anti-parallel to 
Vc, and k has a finite component along Bq. In 
principle, each configuration (vc, k, Bq) requires a 
separate treatment; however, changing the direc¬ 
tion of Bq only inverts the direction of propaga¬ 
tion and does not produce a new dispersion rela¬ 
tion. Therefore, we introduce the dimensionless 
core electron drift 5vc = |uc/r’A|sign(k • Vc), where 
VA = 'JBq!A irnpUip is the proton Alfven speed. 
The values of the remaining parameters are given 
in Tabled! 

We study the dispersion relations of KAWs 
by numerically solving the full set of the lin¬ 
ear Maxwell-Vlasov equations. Linear Maxwell- 



Fig. 1.— Schematic diagram for the plasma con¬ 
sidered in this letter. We annotate the parallel 
wave phase speed Vres = Wr-/fc[| for our resonance 
analysis in Section |31 
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Vlasov theory has been desc ribed in detail i n the 
liter a ture ( see, for instance, SwansorJ ( IQSflh and 
Stixl (1993)). We summarize the relevant results 


here. Linear plasma waves are eigenmodes of the 
wave equation: 


k X (k X E) + • E = 0, (1) 

d- 


where E is the fluctuating electric field in Fourier 
space. The dielectric tensor e incorporates con¬ 
tributions from each plasma component: e = 
1 -I- EsXj,, where Xs is the susceptibility ten¬ 
sor for species s. We obtain from the lin¬ 
earized and Fourier-transformed Vlasov equation 
and Maxwell’s equations. 


To study the contribution to wave damping 
(growth) from each component, we calcula te the 
relative particle heating rates as ( Stix Il993) 


P, = E*.x:|„-E/41F. (2) 

where x1 = (Xs — x|) /2f is the anti-Hermitian 
part of Xs, Wr = 3? (w) and W is the wave energy. 

In the weak damping limit, Pg gives the frac¬ 
tion of wave energy damped by plasma compo¬ 
nent s during one wave period. The total heating 
rate Ptotai = corresponds to the wave en¬ 

ergy dissipation rate directly calculated from the 
damping rate, i.e., 1 — where uji = 3 (w) 

and T = 27: jujr is the wave period, as long as uji is 
small compared to Ur- In this letter, we focus on 
the linear damping of KAWs. Electron drifts may, 
however, alter wave-particle interactions to such 
a degree that the corresponding waves become un¬ 
stable, which we note for the sake of completeness. 
Notice that Pg may become negative, correspond¬ 
ing to cooling of species s. 


Quantity 

Value 

Proton plasma beta /3p 

0.4 

Temperature ratio Tc/Tp 

2 

ThiTp 

10 

Density ratio Uolnh 

9 

Proton bulk drift VpjvA 

0 

Angle between k, Bq 

89° or 91° 


Table 1: Values of fixed plasma parameters in our 
study. 


3. Results 


Figures [IKa)-(b) show dispersion relations of 
KAWs in plasmas with 5vc G [—4,4]. Panel (a) 
demonstrates that electron drifts have little effect 
on the re al part of the KAW frequencies, consis¬ 
tent with Gary et al. ( 1975[) . Indeed, the protons 
rather than the electrons mainly determine the dy¬ 
namics of low-frequency KAWs. Panel (b) com¬ 
pares damping rates of KAWs in plasmas with dif¬ 
ferent electron drifts. Negative (positive) 5vc lead 
to an enhancement (reduction) of the KAW damp¬ 
ing at all wavenumbers. Panel (b) shows that 
5vc —4 leads to an increase in the damping rate 
by ~ 50% compared to the case without electron 
drifts. For sufficiently large positive (5fc, KAWs 
become unstable in certain wavenumber ranges. 
This insta bility is an examp le of a “heat flux in¬ 
stability” (jGarv et al.lll975h . The two-population 
electrons in our model introduce an electron heat 
flux (third moment of the electron VDF), which 
may provide energy to drive certain wave modes 
(e.g. . whistler, Alfven , and magnet osonic) unsta¬ 


ble ( Gary et al. 1975 . 1998al Il999l l. These heat 


flux instabilities, in return, regulate the electron 
heat flux and hence the electron drifts. We note 
that there is no conflict between the damping and 
growth of KAWs and other wave modes. For in¬ 
stance, parallel whistler waves, whose instability 
threshold is much lower than that of the Alfvenic 
instability, are unstable even if 5vc > —4. There¬ 
fore, the Alfvenic instability cannot regulate the 
electron heat flux in the solar wind and is more of 
academic interest under typical solar-wind condi¬ 
tions since the instability with the lowest threshold 
constrains the electron heat flux once triggered. 
Figure He) compares the total particle heating 
rate (solid) with the wave energy dissipation rate 
(dotted). When 1 — < 0.5, i.e., the wave 

retains more than half of its energy after a wave 
period, Ptotai ~ I — Hence we show that 

Pg is a good measure of energy flow as long as 
damping is weak. 


Figures [3|(a)-(d) present heating rates for pro¬ 
tons, all electrons (core -I- halo), as well as sep¬ 
arately for core electrons and halo electrons. Re¬ 
call that positive Pg indicates heating for species s, 
and that negative Pg indicates cooling. Panel (a) 
shows that protons are always heated, although 
electron drifts can slightly modify the values of 
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Fig. 2.— (a)-(b): Dispersion relation of KAWs 
in plasmas for different 5vc- Wr and uii are the 
real and imaginary parts of the wave frequency. 
Pp = VTp/^p is the proton gyro-radius, VTp = 
y^2kBTp/mp is the proton thermal speed, and 
flp = eBo/nipC is the proton gyro-frequency, (c): 
Total particle heating rates {Ptotai , solid lines) and 
wave energy dissipation rates (1 — dots) as 

functions of wavenumber. 


the proton heating rate. On the other hand, pan¬ 
els (b)-(d) show that an electron population can 
experience both cooling and heating depending on 
5vc- While both core and halo are heated in the 
absence of bulk drifts, sufficiently large positive 
5vc lead to core cooling and at the same time sig¬ 
nificantly enhance halo heating. Negative Svc have 
the opposite effect on electron heating. In gen¬ 
eral, electron drifts have a larger effect on electron 
heating than on proton heating. Electron drifts 
modify the damping rates of KAWs by changing 
the efficiency of wave-particle interactions for the 
different electron populations. 


Figure |4] shows how the energy partition de¬ 
pends on electron drifts at k±pp = 1. Pc/Pp and 
Ph/Pp vary wildly with 5vc- However, since elec¬ 
tron drifts affect core heating and halo heating 
in opposite ways, the total electron heating and 
consequently the electron-proton energy partition 
reveal a more moderate dependence on electron 
drifts. Nevertheless, a core drift of 5vc ^ —4 in¬ 
creases Pe/Pp by ~ 25%. Without electron drifts, 
electron heating dominates over proton heating in 
hydro gen plasmas at fc_LPp ^ 1 for (3^1 ( Quataeii] 
Il998ll . We see in Figure 0] that positive 5vc sig¬ 
nificantly reduces electron-to-proton heating ra¬ 
tio. At 5vc 3, electron heating becomes neg¬ 
ligible compared to proton heating. We note at 
this point that the dependence of energy partition 
on electron drifts is qualitatively the same at other 
wavenumbers, which can be inferred from Figure 

El 


We interpret the dependence of wave damping 
and particle heating on electron drifts in terms of 
a simplified resonance analysis. Since Ur is almost 



<5vc 


Fig. 4.— Ratio of energy partition between elec¬ 
trons and protons as a function of electron drift. 
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Fig. 3.— Dependence of Pproton (a), Peiectron (b), Pcore{c) and Phaio (d) on k±pp in plasmas with varying 
Svc- 


independent of electron drifts, the proton heating 
rate shows only a small variation. However, core 
and halo electron bulk drifts significantly change 
the value and the gradient of the eVDF at the 
resonance speed. 


OJr — nOe 



(3) 


where Hg is the electron gyro-frequency and n is an 
integer dependi ng on the polar ization properties of 
the wave mode ( Marschl 200dl . 


In our case, the most relevant resonance for 
wave-particle interactions with KAWs is the Lan¬ 
dau resonance with n = 0. Electrons with the 
(field-parallel) speed Vres interact resonantly with 
the corresponding wave mode. We illustrate the 
case for negative 6vc schematically in Figure [1] 


The gradient of the distribution functions at the 
resonance speed determines if the particle species 
gains or loses energy, i.e., the signs of Pp, P^ and 
Ph- In the example shown in Figure [l] the res¬ 
onance occurs where the halo VDF has a posi¬ 
tive gradient, and hence the KAW removes en¬ 
ergy from the halo electrons, leading to a negative 


value for Ph- In contrast, core electrons and pro¬ 
tons have resonances where the gradients of their 
VDFs are negative, leading to positive values for 
Pc and Pp. Similarly, the dependence of Ps on 5vc 
can also be inferred from the positions of the reso¬ 
nance points with respect to the VDFs. However, 
a qualitative resonance analysis cannot provide a 
thorough comparison across the plasma compo¬ 
nents. We Hnd that Ps^ defined in Eq. gives 
a better understanding of the energy transfer be¬ 
tween waves and different plasma components. 


It is worth noting that electron drifts seem to 
provide a collisionless energy coupling between 
core and halo electrons (see the case with 5vc = — 2 
for instance). Figure [2Kc) shows that I — « 

Pc + Ph 


Pp. 


With Pc,Pp > 0 and Ph < 0 


(from Figure E]), (l - -I- {-Ph) ~ Pc + T’p 

, meaning that energy flows from damped waves 
and halo electrons into core electrons and protons. 
A similar energy transp ort among electron c ompo- 
nents was suggested bv iGarv et al. (1998b). Ac¬ 
cording to this study, the Alfvenic heat-flux insta¬ 
bility transfers energy from the drifting halo into 
the electron core through the Landau-resonant ab- 
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sorption of the waves, leading to a lower limit for 
/3c. 


4. Discussion and conclusion 


For the sake of clarity, we present results for a 
single propagation angle in this letter. A care¬ 
ful analysis shows that other oblique angles do 
not lead to any qualitative differences. Given a 
fixed KAW power spectrum, we anticipate that 
the presence of electron drifts changes the wave- 
energy dissipation rate, particle heating rates, and 
the proton-electron energy partition ratio signifi¬ 
cantly. 

In the context of the solar wind, since core 
electrons are nearly always observed to drift sun¬ 
ward in the proton frame, positive (negative) 6vc 
correspond to KAWs propagating sunward (anti- 
sunward). Therefore, anti-sunward KAWs expe¬ 
rience stronger damping in the presence of elec¬ 
tron drifts and lead to a stronger core heating, 
halo electron cooling, and a stronger total elec¬ 
tron heating. Since the total energy flux typically 
points anti-sunward in the solar wind, we sug¬ 
gest that electron drifts lead to stronger electron 
heating than expected from previous calculations 
( QuataertlllQQSt iLeamon et al. 1999ll . An inter¬ 
esting corollary of this work is that, since the core 
electron drift spe ed depends on the collisionality 
of the solar wind ( Puluoa et al. 2014h . the heating 
rate as calculated here may therfore also indirectly 
depend on solar wind Couloumb collisions. 

Our model includes relative drifts in the elec¬ 
tron VDF, while it ignores several additional fea¬ 
tures of solar wind electrons, namely the pres¬ 
ence of superthermal electrons, one-sided beams 
(strahl) and temperature anisotropies. By using 
a Maxwellian halo, we ignore high energy tails in 
the electron VDF. However, this assumption does 
not significantly affect the gradients of the electron 
distribution function at v^es ^ va which mainly 
determine the damping. We can apply a similar 
argument to strahl electrons, which also occupy 
a different region in the veloc ity space. Regard ¬ 


ing temperature anisotropies, iGarv et al.l (1197511 


studied the dependence of the core drift threshold 
for the Alfvenic instability on electron tempera¬ 
ture anisotropy and found that the dependence is 
weak in plasmas with f3p < 0.25. We expect a 
similar weak relation in our case, since Landau- 


resonant instabil ities show such a weak depen - 
dence in general (jVerscharen fc Chandran 11201,ih . 
since the strength of this type of resonant inter¬ 
action is determined by the parallel gradients of 
the electron distribution function. If the parallel 
temperature is kept constant, the introduction of 
temperature anisotropy only changes the perpen¬ 
dicular gradients of the VDF and therefore does 
not significantly alter the wave-particle interac¬ 
tion. 


In this letter, we choose a fiducial set of repre¬ 
sentative solar wind plasma parameters to demon¬ 
strate the effects of electron drifts on the damping 
of KAWs. In order to fully account for all of the 
relevant effects in the solar wind, it is important 
to conduct a full scan of the corresponding param¬ 
eter space. This endeavor is beyond the scope of 
this work and will be presented in a future paper. 

Our work shows that, despite our limiting as¬ 
sumption of a superposition of Maxwellians to rep¬ 
resent the electron distribution function, its fine- 
structure has a strong influence on the propaga¬ 
tion and damping properties of KAWs. Therefore, 
our work is of relevance in the broader context 
of all collisionless astrophysical plasmas in which 
non-Maxwellian electron distributions can develop 
and persist. Considering that most astrophysical 
plasmas are in a turbulent state and that KAW- 
turbulence is believed to be the dominant type of 
plasma turbulence on small scales, our work sug¬ 
gests that the effects of non-Maxwellian electron 
distributions be carefully accounted for in studies 
of collisionless astrophysical plasmas in general. 
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NASA grant APL-975268. C. H. K. Chen is sup¬ 
ported by an Imperial College Junior Research Fel¬ 
lowship. C. Salem is supported by NASA grant 
NNXI4AC07G. D. Verscharen is supported by 
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